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Abstract. Resonant vibrational and rotation-vibration excitation cross sections for 
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exchange plus polarization calculations performed using the R-matrix method are used 
to estimate resonance positions and widths as functions of intcrnuclear separation. 
The effects of nuclear motion are considered using a local complex potential model. 
Good agreement is obtained with available experimental data on excitation from the 
vibrational ground state. Excitation rates and cross sections are provided as a functions 
of the initial CO vibrational state for all ground state vibrational levels. 
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1. Introduction 

When a spacecraft enters an atmosphere at speeds exceeding the local speed of sound, 
a shock wave is formed behind it and its kinetic energy is transformed into heat. The 
energy delivered to the gas in this process promotes excitation of the molecular internal 
degrees of freedom (rotational, vibrational and electronic) and chemical reactions 
(including dissociation and ionization). The hot reacting (and radiating) gas is often 
in thermal and chemical nonequilibrium conditions and it is this complex system that 
interacts with the vehicle surface. 

In this paper we consider electron-impact resonant vibrational excitation and 
rotational- vibrational excitation of carbon monoxide over the entire range of vibrational 
levels supported by its ground electronic state, CO(X 1 S + ), and over an extended 
electron temperature range. The energy range of concern here is 0-10 eV. This data 
provides important input information into spacecraft entry models for atmospheres such 
as those on Mars and Venus as well as to study other processes in the atmospheres 
of these planets [1] and in comets [2]. It is also useful for understanding processes 
involved in the CO laser [HI H] and to study CO plasma in presence of electrical 
discharge E]. Infrared emission from CO in the upper atmospheres of Mars, Venus 
and several other planets is a subject of current theoretical and experimental interest. 
The first measurements of cross sections for low-energy electron impact excitation of 
the vibrational levels of the ground state of CO have been made by Schulz [7J. Recent 
new measurements showing the contribution of electron impact relative to emissions by 
other mechanisms, have been reported pQ. 

Direct vibrational excitation of diatomic molecules by electron impact is, in general, 
an inefficient process because of the small electron-to-molecule mass ratio. However, 
when the incident electron can attach to form a temporary negative ion, vibrational 
excitation cross section can be hugely enhanced [8]. These processes are called resonant 
collisions. In this paper the following reaction is treated: 

e" + CO(V 1 S+,t; i ) -> CCr( 2 n) -> e'~ + CO(X l Z + ,v f ) , (1) 

where Vi and Vf are the initial and final vibrational levels of the molecule. A number of 
reviews on this subject are available [9j []JJ [UJ [T2] . 

Resonance enhanced vibrational excitation of CO from its vibrational ground 
state (vi = 0) has been well-studied experimentally, notably in recent work by 
Allan [13], Poparic et al. [H] and by Gibson et al. [15]; earlier experimental work is 
reviewed by Brunger and Buckman [12]. Theoretically the best calculations are due 
to Morgan [T6t ITT] who used an R-matrix method to characterize the resonance as a 
function of CO internuclear distance. As will be shown below, her studies show rather 
good agreement with both data available to her at the time and subsequent studies. 
The various studies show that electron collisions involving the well-known, low-lying 2 I1 
shape resonance lead not only to a large vibrational excitation cross section but also 
show that Av (= Vf — v,j) can be large: excitation cross sections with Av = 11 were 
measured by Allan [T3] . 
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So far studies of electron impact vibrational excitation of CO have concentrated 
almost exclusively on excitations from the vibrational ground state. However for 
modeling CO in hot environments it is necessary to know the corresponding cross 
sections starting from excited vibrational states, Vi > 0. In this work we develop a 
model which reproduces the known data for V{ = and then use it to obtain vibrational 
excitation cross sections and rates for the whole range of possible initial vibrational 
states. 

The paper is organized as follows: In section |2] the theoretical model is presented 
for both electron dynamics and nuclear motion; in section [3] the calculated results are 
presented and in section H] they are discussed. 

2. Theoretical model 

2.1. Electron collisions 

The R-matrix method |18l [19] is used to obtain a complex potential energy curve for 
the resonance. Our methodology and calculations are heavily based on the work of 
Morgan [16] who, as is shown below, obtained excellent results for the excitation process. 
For the present study it was necessary to extend the range of Morgan's calculations to 
allow for excitation to high-lying bound states. To do this we follow the prescription 
for the electron collision calculations given in her paper but, since the paper does not 
provide entire details of the calculation, it was necessary to first perform a number of 
test studies. The nuclear dynamics is treated in the local-complex-potential model as 
explained in the next subsection. 

Calculations were performed using the diatomic molecule implementation of the 
UK Molecular R-matrix codes [20] which uses Slater Type Orbitals (STOs) to represent 
the target and numerical functions to represent the continuum. We used the STO 
basis given by Kirby-Docken and Liu [2T] to generate a total of 16a and 12ir molecular 
orbitals, and numerical Bessel functions for partial waves up to / = 6 for the continuum 
inside the R-matrix sphere which had a radius of 10 ao- To avoid linear dependence 
problems, two of the continuum 7r orbitals were removed by Lagrange orthogonalization 
to the target orbitals [22]. Resonance positions and widths were obtained by fitting the 
2 I1 eigenphase sum at geometry to a Breit-Wigner form using program RESON [23J. 
The resonance becomes very broad at short internuclear separations and, under these 
circumstances, the background eigenphase can vary significantly across the resonance 
resulting in fits that are less stable. The effect of this can be seen in the behavior of the 
fitted width as a function of CO bondlength. See Fig. [2] below. 

Morgan tested a number of models but found that a Static Exchange plus 
Polarisation (SEP) model performed best. This model uses a Hartree-Fock (HF) target 
wavefunction and includes polarisation effects augmenting the scattering wavefunction 
with so called two particle - one hole (2p-lh) configurations. These configurations 
involve simultaneously exciting a single target electron into an unoccupied target 
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("virtual") orbital and placing the scattering electron in a virtual orbital. Results for 
resonances parameters are well-known to be sensitive to the precise choice of parameters 
in an SEP calculation [21]. Our final model froze the C and O Is and 2s electrons 
and considered excitation of remain 6 "valence" electrons into virtual orbitals, all of 
which were retained. At the CO equilibrium bondlength this model gives a resonance 
position and width of (1.67, 0.82) eV which can be compared to Morgan's values of (1.68, 
0.95) eV. It can be seen that the positions are in excellent agreement but our width is 
somewhat narrower. As discussed below, it was decided to increase our calculated widths 
by 10 % to improve agreement with experiment. 



2.2. Nuclear motion 

Our treatment of vibrational excitation by low-energy electron impact follows the 
general formulation given in the refs. [H [25]. At low energies the cross section is 
dominated by negative ion resonance contributions, that is the incident electron is 
temporarily captured by the molecule and a negative intermediate state occurs, then 
the resonant state decays into a new state. In this paper, process (OQ) is described 
within the framework of the local- complex-potential (LCP) model [261 123 I2H], that is 
an approximation to the more general non-local theory, which is appropiate when the 
resonance width is much larger than the spacing between the target vibrational levels (for 
CO ~ 0.1 eV). In general the resonance width depends on the energy as well as on the 
internuclear distance, T(E, R); the LCP approximation replaces this energy-dependence 
with the value E(R), the fixed-nuclei resonance energy [29] : 

T(E,R)^T(E(R),R) = T(R). (2) 

In the following all dynamical quantities are understood energy- independent. 

In the LCP approach the electron-molecule cross section, in the rest frame of the 
molecule for vibrational transition Vi — > Vf and for an incident electron with energy 
is given by [8]: 

^ = Kx/|v|£>l 2 , (3) 

where m e is the electron mass and g = 2 is a statistical factor, k^,^ = 2 m e / h 2 is the 
incoming (outgoing) electron momenta, £(R), depending on the internuclear distance 
R, is the solution of the nuclear wave equation of the negative ion, with total energy 
E — €i + e Vi , 

Tj + V'+A-'-r-E^J Z(R) = -V Xi(R) , (4) 

and Xi(f)(R) is t ne initial (final) vibrational eigenfunctions of the neutral molecule with 
eigenvalue e v . (f) , given by the equation: 

(Tj + V )xn(R) = e VnX n(R). (5) 
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The kinetic energy operator, Tj(R), including the centrifugal potential, is expressed as 

T (R) - * ^ + J ( J+1 ^ m 

where \x is the nuclei reduced mass and J the angular momentum quantum number 
of the target molecule. The adiabatic potentials of the target and resonant electronic 
states, Vo(R) and V~(R), have been represented by a Morse function: 

V(R) = D e [l- e - a{R - R ^] 2 + W . (7) 

In the LCP model, the discrete-state-continuum coupling potential V is expressed 

as: 

where k(R) is defined by 

#{R) = ^\V--V Q \. (9) 

The level shift operator A(R) and the resonance width T(R) are discussed in the next 
section. 

Once the collision cross section <7j/ and the electron energy distribution are known, 
the vibrational excitation rate coefficient Kjj can be evaluate. Assuming that the 
electron energy distribution is Maxwellian, the rate coefficient can be expressed, as 
a function of the electron temperature T, as: 



9 f°° 
K if (T) = -=( K T)- 1 - 6 / ea lf (e)e-^ T de 



(10) 



where k is the Boltzmann constant and eth = Vf — Vi is the threshold energy for the 
process 



3. Results 

The ground state of the potential energy curve for carbon monoxide CO(X was 
calculated using Molpro [50] and an aug-cc-pwCV5Z GTO basis in a Davidson-corrected 
Multi-Reference Configuration Interaction (MRCI) calculation. The calculated points 
have been fitted with a Morse function and the parameters are displayed in Table HJ 
compared with those of different authors. The CO~( 2 II) potential curve was taken 
from the R-Matrix described above up to the crossing point with the CO potential at 
R = 2.6 ao; at longer bondlengths a Molpro calculation was used to extend the curve 
to its asymptotic limit. This curve matched the resonance curve at the crossing point. 
The resonance potential was made to reflect the true shape of the CO potential energy 
curve by using resonances energies relative to our calculated target curve rather than 
the HF curve used in the calculations. Our curves give an asymptotic electron affinity 
of 1.43 eV, close to the observed oxygen atom electron affinity of 1.46 eV [31] . It is 
important to get these details correct as the vibrational excitation cross sections are 
very sensitive to the relative positions of the CO and CO - curves. Fig. [T] shows in 



Resonant VE electron-CO cross sections 



6 





This work 


CO(X 1 £+) 
Theory [32] 


Exp. [33] 


CCT 
This work 


( 2 IT) 

Theory [16] 




(eV) 


11.26 


11.19 


11.22 


9.76 




Re ( 


a.u.) 


2.13 


2.13 


2.13 


2.30 


2.25 


w 


(eV) 









1.49 





Table 1. Dissociation energy {D e ), equilibrium distance (R e ) and minimum position 
(W) for the CO and CO~ Morse-like parameters of Eq. (J7J, compared with those of 
different authors. 




Internuclear distance (a.u.) Internuclear distance (a.u.) 

Figure 1. Left: CO and CO~ potential energy curves. Right: CO ground state 
potential energy curves including the rotational contribution coming from different 
values of the rotational angular momentum, </, as indicated in the figure. 

the left panel the obtained CO and CO~ potential energy curves as a function of the 
internuclear distance. The CO(X 1 E + ) ground state has been found to support 81 
vibrational states. The righthand panel of Fig. 1 shows curves representing the ground 
state potential energies plus the angular contribution coming from different rotational 
states (J = 0, J = 100, and J = 200). 

The resonance width, T(R), has been fitted with a polynomial function, 

T(R) = (ci + c 2 R + c 3 R 2 + c 4 i? 3 + c 5 R*) 9 [V~{R) - V (R)] , (11) 

where c 1 = -302.66 eV, c 2 = 635.8 eV/a , c 3 = -480.06 eV/ag, c 4 = 156.9 eV/ag, 
c 5 = —18.88 eV/a^; 9 is the step function. This fit and the calculated results are 
displayed in Fig. [2j In the local version of complex-potential model is not possible to 
calculate the level shift A(i?) from first principles (this requires T(E,R)). In the LCP 
approach the level shift is therefore treated a phenomenological parameter external to 
the model. We use this degree of freedom to match correctly the position of the peaks 
using the 4 1 and — > 2 cross sections of Allan [13] . Moreover we scale the width 
T(i?) to correctly reproduce the experimental height of the main peaks. In practice 
T(i?) was multiplied by 1.1 and the A parameter was set to 0.035 eV. Fig. |3] compares 
our ab initio cross sections and the final results. 

Fig. [3] compares our vibrational excitation cross sections for transition starting 
from the lowest vibrational level (vi = 0), with the measurements of Allan [13] and 
the previous ab initio predictions of Morgan [16]. It can be seen that the agreement 
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Figure 2. Resonance width, T(R), as a function of bondlength. The figure shows the 
the R- matrix points and the fitting curve given by Eq. 
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Figure 3. e-CO resonant cross sections compared with the experimental results of 
Allan [13] and the theoretical R- matrix calculation of Morgan [16] . 



between the three studies is very good for the lower vibrational levels. The discrepancy 
for high-lying levels, as shown by — > 10 even when the experimental error of 20% is 
taken into account, suggests a full non-local model and an energy-dependant widths, T, 
are needed for these states. However, we note that the cross sections for these extreme 
excitations are very small. 
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Figure 4. c-CO resonant cross sections as a function of the incident electron energy. 
Upper panels: elastic processes involving the u, — v f vibrational levels, as indicated 
in the figures. Lower panels: Vi — > Vf processes starting from Vi = 1 (left), Vi = 10 
(right) and vt > u,. 



Having developed a satisfactory model for electron impact vibrational excitation of 
CO we apply it to the whole range of vibrational states supported by the CO molecule. 
Figure H] gives sample results showing that the excitation cross sections are strongly 
state dependent. 

Fig. [5] shows the rate coefficients for the inelastic transitions Vi — > Vf, calculated 
from Eq. (fTOj) assuming a Maxwellian electron energy distribution function. The left 
panel of this figure illustrates the behavior of the rate coefficients for the inelastic 
processes starting from the level Vi — 0. The decrease of the rates with Vf can be 
attributed to the reduction in the corresponding cross sections, as is seen in Fig. |3j 
Analogous behavior is observed in the right panel, where the rates for vibrational 
excitations starting from Vi = 10 are shown. 

We have also tested the effect of including rotational motion in the calculation. 
Figs. M and [7] show, respectively, some examples of cross sections and rates coefficients 
for different values of the angular momentum quantum number J. For cool CO samples, 
which only probe low-lying rotational states, the effect of including rotational motion 
is small. However for hot samples of CO, which occupy highly excited J states, the 
magnitude, structure and position of the resonance enhanced excitation cross section 
changes significantly with rotational excitation. 
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Electron temperature (eV) Electron temperature (eV) 

Figure 5. e-CO resonant vibrational-inelastic excitation rate coefficients for the 
processes — > Vf (left panel) and 10 — > Vf (right panel). 
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Figure 6. e-CO resonant vibrational excitation cross sections calculated for different 
values of J. 

4. Summary 

The aim of this work is to provide electron impact excitation cross sections and rate 
coefficients of CO for modeling purposes. The theoretical cross section calculations 
have been performed in the framework of the local complex potential model for resonant 
collisions while the input parameters, adiabatic potential energies and widths, have been 
computed by the R-matrix method. The results obtained are in good agreement with 
the existing experimental and theoretical data. 

A full set of cross sections and state-dependent rates can be obtained from [34] . 
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Figure 7. Rate coefficients as a function of the eiectron temperature for e-CO elastic 
scattering calculated for J ~ 0, 100 and 200. 
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